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Abstract: This paper studies the nonlinear stochastic partial differential equation 
of fractional orders both in space and time variables: 






p(u{t,x))W(t,x) , 


t > 0, X e 


where W is the space-time white noise, a € (0,2], f3 € (0,2), 7 > 0 and ^ > 0. 
Fundamental solutions and their properties, in particular the nonnegativity, are 
derived. The existence and uniqueness of solution together with the moment bounds 
of the solution are obtained under Dalang’s condition: d < 2a min (27 — 1, 0). In 
some cases, the initial data can be measures. When f3 e (0,1], we prove the sample 
path regularity of the solution. 
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1 Introduction 


In this paper, we will study the following nonlinear stochastic time-fractional diffusion equa¬ 
tions: 


< m( 0, ■) = h 

d 

M(0,-)=ho, ^M(0,-)=/il 


p{u{t,x))W{t,x) , 


t > 0, X e 
if/3e ( 1 / 2 , 1 ], 

if/Se (1,2), 


( 1 . 1 ) 


with a G (0, 2] and 7 > 0. In this equation, A = is the Laplacian with respect 

to the space variables and (—A)"/^ is the fractional Laplacian. W denotes the space-time 
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white noise, z/ > 0 is the diffusion parameter. The initial data /r, /xq and /ii are assumed to 
be some measures, p is a Lipschitz continuous function, denotes the Caputo fractional 
differential operator: 




dP{t) := 


_ I r(m - 13 ) (i - r)<s+‘ 


di’ 


:/(i) 


ifm — l</3<m, 
if /S = m , 


and If is the Riemann-Liouville fractional integral of order 7 > 0: 


ifm ■■= 


r(7) 


{t — sy ^/(s)ds, for f > 0 , 


with the convention If = Id (the identity operator). We refer to [16, 29, 30] for more details 
of these fractional differential operators. 

This paper is an extension of a recent work by the first author [4], where the case 7 = 
\ff\ — y a = 2 and d = 1 is studied. Here, [/S] is the smallest integer not less than The 
interested reader can find motivations of the model in that reference. The fractional integral 
operator If smooths the noise term. Removing this integral operator by setting 7 = 0, one 
may expect that the solution becomes less regular. As proved in [4], when a = 2, d = 1 and 
7 = Iff] — /3, there is a mild solution for all ft G (0, 2). This is no longer true if this fractional 
integral operator is not there. In particular, we will show that, when a = 2, d = 1 and 
7 = 0, the mild solution exists only for f G (2/3,2) instead, which is a direct consequence 
of condition (1.9) below. 

Motivations for stochastic partial differential equations (spde) with time-fractional deriva¬ 
tive can be found in many recent papers [4, 11, 18, 26]. For convenience, we call equation 
(1.1) with f G (0,1] the slow diffusion equation, and equation (1.1) with f G (1,2) the fast 
diffusion equation. 

When d = l, /3 = 2, a = 2 and 7 = 0 , the spde (1.1) reduces to the stochastic wave 
equation (SWE) on M: 


dff 


V 




2 dx'^ 


u{t,x) = p{u{t,x))W{t,x) 


( 1 . 2 ) 


with the speed of wave propagation { 0 / 2 )^/“^. When d=l, fi = l, a = 2 and 7 = 0, the 
spde (1.1) reduces to the stochastic heat equation (SHE) on M: 




V 




2 dx‘^ 


u{t,x) = p{u{t,x))W{t,x) 


(1.3) 


These two special cases have been studied carefully; see [3, 5, 6 , 7, 12]. The spde (1.1) for 
ft G (0,1] and u = 1 — ft has been recently studied in [25, 26]. When the noise does not 
depend on time, a similar model with a general elliptic operator has been studied in [18]. 
Another related equation is the stochastic fractional heat equation (SFHE) on M: 

- xDf^ u{t,x) = p{u{t,x))W{t,x) , 
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(1.4) 








which has been studied recently in [8, 10]; see also [15, 19]. 

All investigations on spde’s of the above kinds require a good study of the corresponding 
Green functions. As proved below, there is a triplet 

{ Z{t, x), Z*{t, x) , Y{t, x) : (t, x) G M+ x } , 

depending on the parameters {a, /3, 7 , u), such that the solution to (1.1) with p{u{t, x))W{t, x) 
replaced by a nice function f{t,x) is represented by 

(Z(t, ■) */i)(a;) + (F x/) (t,x), if/? G (0,1], 

u{t,x)={ (1.5) 

{Z*{t,-) * fio){x) + {Z{t,-) * pi){x) + (Fx/) {t,x), if /? G (1,2), 

where denotes the convolution in the space variable: 


(Z(t,-)tii)(x):= Z{t,x-y)li(dy), 


( 1 . 6 ) 


and “x” denotes the convolution in both space and time variables: 


{Y-kf){t,x):= / Y{t-s,x-y)f{s,y)dsdy. 

Jo Jm.<i 

These fundamental solutions are expressed using the Fox H-function [20]. 

If we denote the solution to the homogeneous equation of (1.1) by Jo{t,x), i.e., 

^ f (^(^1 •) */^) (a^) if/?G(0,l], 

\{Z*{t,-)*po){x) + {Z{t,-)*pi){x) if/? G (1,2), 

then the rigorous meaning of (1.1) is the following stochastic integral equation: 
u(t,x) = Jo(t,x) + I (t,x), where 


Iit,x)= I Y{t- s,x-y)p{u{s,y))W{ds,dy). 

'[0,t]xR 


(1.7) 


( 1 . 8 ) 


The stochastic integral in the above equation is in the sense of Walsh 

To establish the the existence and uniqueness of random field solutions to (1.1), the first 
step is to check Dalang’s condition [14]: 

[ ds [ dy |F(s, ?/)p < cx), for alH > 0, 


which is equivalent to the following condition (see Lemma 5.3): 


a 


d < 2a + — min (27 — 1, 0) =: 0, 

y 


(1.9) 
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which is equivalent to 


/;+7>tii + ^ 

a 


and d < 2a. 


( 1 . 10 ) 


Note that (1.10) implies that the space dimension should be less than or equal to 3. Among 
all possible cases in (1.9), the following two special cases have better properties: 


7 = 0 or a > d = 1, 
a > d = 1, 


( 1 . 11 ) 

( 1 . 12 ) 


As shown in Lemma 4.3 and Remark 4.4 below, under both conditions (1.9) and (1.11), the 
function Y (1, x) is bounded at a; = 0. Moreover, under (1.9) and (1.12), all functions Z{l,x), 
Z*(l,x) and F(l,x) are bounded at a; = 0. 

We prove the existence and uniqueness of random held solutions to (1.8) in the following 
three cases: 

Case I: If we assume only Dalang’s condition (1.9), we prove the existence and uniqueness 
when the initial data are such that 


sup |Jo(-S; 2 ;)| < oo, for all f > 0, 

(s,x)£[0,t]xR‘^ 


(1.13) 


which is satished, for example, when initial data are bounded measurable functions. 

Case II: Under both (1.9) and (1.11), we obtain moment formulas that are similar to those 
in [7, 4, 8]. The initial data satisfy (1.13). 

Case III: Under both (1.9) and (1.12), the initial data can be measures. Let AI(M) be the 
set of signed (regular) Borel measures on M. For a; G M, dehne an auxiliary function 


'■= exp (-7 . 


(1.14) 


where L/3J is the largest integer not greater than /?. Note that the difference between [/S] and 
[/9J + 1 for /3 G (0,2) is only at /d = 1. The initial data are assumed to be Borel measures 
such that 


(IfI * < oo, 

f 1 

sup / \fi\{dx) 
ysM Jr 


1 + |x — 


for all 7 > 0 and a; G M, if a = 2, 
< +CX), if a G (1, 2), 


(1.16) 


where for any Borel measure /a, /a = /a+ — /a_ is the the Jordan decomposition and |/a| = 
/a+ + /a_. We use A4a,/3(M) to denote these measures. In this case, we prove the existence 
and uniqueness of a solution to (1.4) for all initial data from AJa,/ 3 (lR). 

Here are some special cases: 

(1) For (1.3), i.e., a = 2, (3 = 1 and 7 = 0, the set of admissible initial data studied in [7] is 
A4 h(M), which corresponds to AJ 2 ,i(l^) in this paper. 
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(2) Under the condition that d = 1, a = 2, /3 G (0, 2), 7 = \(3'\— (3 (as in [4]), one can easily 
verify that condition (1.9) is always true. The possible initial data is A4^(M), which 
corresponds to Al 2 ,/ 3 ( 1 ^) in this paper. 

(3) If 7 = 1 — /3 and (3 G (0,1), then it is ready to see that (1.9) reduces to 

d < Q!min(2,/3“^), 

which recovers the condition by Mijena and Nane [26]. 

(4) If 7 = 0, then (1.9) becomes 


Moreover, if a = 2 and d = 1, then this condition becomes (3 > 2/3, which coincides to 
the condition in [11, Section 5.2]. 

(5) If /3 = 1 and 7 = 0, then Dalang’s condition (1.9) reduces to a > d. Since a G (0, 2], we 
have that a G (1, 2] and d = 1, which recovers the condition in [ 8 ]. 

As in [ 6 , 7, 8 ], we will obtain similar moment formulas expressed using a kernel function 
IC(t,x) when (1.11) is satisfied. For the SHE and the SWE, this kernel function IC(t,x) has 
explicit forms. But for the SFHE [ 8 ], (1.1) with < 7 = 1,7 = o; = 2 in [4], and the 

current spde (1.1), we obtain some estimates on it. In particular, we will obtain both upper 
and lower bounds on IC(t,x). 


After establishing the existence and uniqueness of the solution, we will study the sample- 
path regularity for the slow diffusion equations (i.e., the case when (3 G (0,1]). Given a 
subset D C M_,_ x and positive constants /3i, [32, denote by the set of functions 

V : M+ X I— >■ M with the following property: for each compact set K (Z D, there is a finite 
constant C such that for all {t, x) and (s, y) G K, 

\v{t, x) — v{s, y)\ < C [\t — + \x — 2 / 1 ^^) . 


Denote 


^/3i-,/ 32-(-^) ■ (0,/3i) Gcjg (0,02) ^ai,a2{^) ■ 


We will show that for slow diffusion equations, if the initial data has a bounded density, i.e., 
y{dx) = f{x)dx with / G L“(M'^), then 


where 0 is defined in (1.9), and := (0,oo). 


When the initial data are spatially homogeneous (i.e., the initial data are constants), so 
is the solution u(t,x), and then the Lyapunov exponents 


Trip := 


lim sup - logE [|M(f, x) 1 ^] , 

t—KX) t 


(1.17) 
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(1.18) 


m ■= lira inf - log E [\u{t, x) |^], 

^ t—^OO t 

do not depend on the spatial variable. In this case, a solution is called fully intermittent if 
mi = 0 and m 2 > 0 (^ 6 ® Dehnition III.1.1, on p. 55]). As for the weak intermittency, 
there are various dehnitions. For convenience of stating our results, we will call the solution 
weakly intermittent of type / if m 2 > 0; weakly intermittent of type //if m 2 > 0. Clearly, 
the weak intermittency of type I is stronger than the the weak intermittency of type II, but 
weaker than the full intermittency by missing mi = 0. The weak intermittency of type II is 
used in [19]. 

The full intermittency for the SHE and the SFHE are established in [1] and [10], respec¬ 
tively. The weak intermittency of type I and II for SWE are proved in [ 6 ] and [13, Theorem 
2.3], respectively. We will establish the weak intermittency of type II for both slow and fast 
diffusion equations. For some slow diffusion equations, we will prove the weak intermittency 
of type I. Moreover, we show that 

IFp < C . (1-19) 

It reduces to the following special cases: 

(1) The SHE case, i.e., /3 = 1, 0 = 2, 7 = 0 and d = 1: Trip < C p^. See [1, 7, 19]; 

(2) The SWE case, i.e., /3 = 2, 0 = 2, 7 = 0 and d = 1: Trip < C See [ 6 ]; 

(3) The SFHE case, i.e., /3 = 1, 7 = 0 and d = 1: Trip < C p2+i/("-i). See [ 8 ]; 

(4) The time-fractional diffusion equation case as in [4] that a = 2, 7 = [/?]—/3 and d = 1: 

4r/31-/3 

rUp < C p 4 r/ 31 - 2-;3 ■ 

(5) The time-fractional spde as in Theorems 3.11 and 3.12 of [9] with d = 1 and k = 

_ 2a-/3 

urip < C p'^a-p-a when 7 = 0, 

uTLp < C p• 2 <=‘\p^-p-a when 7 = [/d] — (3. 


In order to obtain some lower bounds for the moments, we prove that under the following 
cases 


"Case I: a G (0, 2], //G (0,1), d G N, 7 > 0, 

Case H; a G (0, 2], /? = 1, d G N, 7 G {0} U (1, cx)). 

Case HI: l</d<Q;<2, d<3, 7 > 0 , 

Case IV: 1 < (3 = a < 2, d < 3, 7 > {d + 3)/2 — f3, 


( 1 . 20 ) 


^In [9], the constant k is the exponent for the Riesz kernel, and the case k = 1 corresponds to the 
space-time white noise. 
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the fundamental solution Y is nonnegative (see Theorem 4.6 blow). These results generalize 
those obtained by Mainardi et al [24], Pskhu [28], and Chen et al [9]. In the end, we derive 
some lower bounds for the moments of the solution to (1.4) under (1.9) and (1.20). 

This paper is structured as follows. In Section 2 we hrst give some notation and prelim¬ 
inaries. The main results are stated in Section 3. The fundamental solutions are studied in 
Section 4. The proof of the two existence and uniqueness theorems are given in Section 5. 
Finally, in Appendix, we prove some properties of the Fox H-functions. 

2 Some preliminaries and notation 

Let W = {Wt{A) :AeBb O} be a space-time white noise dehned on a complete 

probability space (fl, J^, P), where Bb (M'^) is the collection of Borel sets with hnite Lebesgue 
measure. Let 

= (T (Ws{A) : 0 < s < t, A G Bb V TV, t > 0, 

be the natural hltration augmented by the cr-held TV generated by all P-null sets in P. We 
use ll'llp to denote the LP(f2)-norm {p > 1). In this setup, W becomes a worthy martin¬ 
gale measure in the sense of Walsh [33], and J/jg X(s, j/)hF(ds, dy) is well-dehned for a 
suitable class of random helds {X(s,j/), {s,y) G M+ x 

Recall that the rigorous meaning of the spde (1.1) is in the integral form (1.8). 

Definition 2.1. A process u = {t,x) G x is called a random field solution 

to (1.1) if 

(1) u is adapted, i.e., for all {t,x) G x u{t,x) is Pi-measurable; 

(2) u is jointly measurable with respect to B {R\ X X P; 

(3) for all {t, x) G x the following space-time convolution is hnite: 

^\\p{u)\\l) {t,x) := [ ds [ d|/V^(f-s,x-|/)||p(M(s ,|/))||2 < Too; 

Jo JRd. 

(4) the function {t,x) h ->■ I{t,x) mapping x into Lf{Vt) is continuous; 

(5) u satishes (1.8) a.s.,for all {t,x) G x R'^. 

Assume that the function p : M i—>■ M is globally Lipschitz continuous with Lipschitz 
constant Lip^ > 0. We need some growth conditions on p assume that for some constants 
Lp> d and ? > 0, 


\p{x)\^ < l?p +x^) , for all X G M. (2.1) 

^This is a consequence of the Lipschitz continuity of p. 
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Sometimes we need a lower bound on p{x)\ assume that for some constants /p > 0 and ? > 0, 

\p{x)f > , for all a; G M . (2.2) 

For all (f, x) G x n G N and A G M, define 
Co (f,x) := Y'^it.x) 

Cn (t, x) := (£o X ■ ■ ■ X Cq) {t, x), for n > 1, {n convolutions), 


/C («, x\ A) := ^ A2<”+‘>£„ (i, A) 


(2.3) 

(2.4) 


n=0 

We will use the following conventions to the kernel functions /C(t,x; A): 

/C(t, x) := /C(t, X] A), /C(t, x) := /C (t, X] Lp) , 

/C(t, x) := 1C (t, X] Ip) , Al!p(t, x) ■.= fC (t, x; Lp), for p > 2 . 

Throughout the paper, denote 

a := 2(1 — /3 — 7 ) + j3d/a. 

Note that 

OL 

(1.9) ^ d<2a + -(2z/-l) ^ 2(/3 + 7 ) - 1 - d/3/a > 0 ^ a < 1. (2.7) 

Let tD°l denote the Riemann-Liouville fractional derivative of order a G M (see, e.g., [29, 
(2.79) and (2.88)]): 


( 2 . 6 ) 


( 2 , 6 ) 


dr 


/M 


iDlfit) := 


r(m —a)df™7o (f — r)“+^" 

dm 


dt"^ 

1 


fit) 


^—a—1 


fis)ds 


Ir(-a) Jo 

We will need the two-parameter Mittag-Leffter function 


\i m — 1 < a <m and a > 0, 

if a = m > 0, (2.8) 

if a < 0. 


Ea,l3iz) 


-r(afc + /?)’ 


a > 0, /3 > 0 , 


(2.9) 


fc =0 


which is a generalization of exponential function, Ei^i{z) = e^; see, e.g., [29, Section 1.2]. 
A function is called completely monotonic if (—l)"'/^’^^(x) > 0 for n = 0,1, 2,...; see [34, 
Dehnition 4]. An important fact [31] concerning the Mittag-Leffler function is that 

X G M+ t—)■ Ea^gi—x) is completely monotonic 0 < a < 1 A /?. (2-10) 

By [20, (2.9.27)], the above Mittag-Leffler function is a special case of the Fox H-function: 

E^^{z) = Hl'l(-z ^ 
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3 Main results 


The first two theorems are about the existence, uniqueness and moment estimates of the 
solutions to (1.1). The second one, in particular, possesses the same form as the one in [4, 
Theorem 3.1]. See also similar results for other equations, e.g., SHE [7, Theorem 2.4], SWE 
[6, Theorem 2.3], and SFHE [8, Theorem 3.1]. 

Theorem 3.1 (Existence, uniqueness and moments (I)). Under (1.9), the spde (1.1) has a 
unique (in the sense of versions) random field solution {u{t,x) : (t,x) G x if the 
initial data are such that 


Ct := sup |Jo(s,x)| < +CX). 

(s,a;)e[0,t] xR'* 


(3.1) 


Moreover, the following statements are true: 


(1) {t,x)^u{t,x) is continuous for all p > 2; 

(2) For all even integers p>2, all t > 0 and x,y ^ 

exp , (3.2) 


\u{t, x) 




2Jl(t,x) + 


7 ' + 2 ^ 


where C is some universal constant not depending on p, and a is defined in (2.6). 


This theorem is proved in Section 5.6. Note that if the initial data are bounded functions, 
then (3.1) is satished. 

Theorem 3.2 (Existence, uniqueness and moments (II)). If Dalang’s condition (1.9) is 
satisfied, then the spde (1.1) has a unique (in the sense of versions) random field solution 
{u{t,x) : {t,x) G X starting from either initial data that satisfy (3.1) under condition 
(1.11) or any BOrel measures from Ada,under condition (1.12). Moreover, the following 
statements are true: 


(1) (t,x)^u(t,x) is U’ifil)-continuous for all p >2; 

(2) For all even integers p >2, all t > 0 and x,y E 




J^{t, x) + ([<^2 +J^] X /C) (t, x), 
2Jl{t,x) + (^[7^+2J2] {t,x), 


ifP = ‘^, 

ifp>2- 


(3.3) 


(3) If p satisfies (2.2), then under (1.9), (1.11) and the first two cases of (1.20), for all 
t > 0 and x,y E W^, it holds that 

\\u{t,x)\\l > Jl{t,x) + ((i^+Jo) x/c) {t,x) . (3.4) 


The proof of this theorem is given in Section 5.5. 

The following theorem gives the Holder continuity of the solution for slow diffusion equa¬ 
tions. 
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Theorem 3.3. Recall that the constants a and 0 are defined in (2.6) and (1.9), respectively. 
If fi ^ (Oj 1] (3.1) holds, then under (1.9), 

sup I |M(t, a;)| Ip <+ 00 , for all T > 0 and p > 2. (3.5) 

(t,3;)e[0,T]xM'^ 

Moreover, we have 

^ C'i(l-(T)-, |min(e-d,2)- ^ (3.6) 


and (1.16) holds. 

Proof. The bound (3.5) is due to (3.1) and (3.3). The proof of (3.6) is straightforward under 
(3.5) and Proposition 5.4 (see [5, Remark 4.6]). □ 

In order to use the moment bounds in (3.3) and (3.4), we need some good estimate on 
the kernel function IC{t,x). Following [4], dehne the following reference kernel functions: 

cy (4z/7rf^) "^^^exp , if a = 2, 


Qa,y{t, x') . < 


Cd td/^ 


^(^2/3/a + |^|2)('='+1)/2’ 


(3.7) 

if a e (0,2), 


for fi E (0, 2) where |a;p = x^-l - \-x‘^, C /3 = 1 if /? e [1, 2) and = 2 R‘^T{d/2)/T{d) 

if /3 G (0,1), and Cd = + l)/2). Dehne also 


Qa,is(t, x'j 



Cd tdl 




if a = 2, 
if a e (0,2). 


(3.8) 


These reference kernels are nonnegative and the constants and q are chosen such that 
the integration of these kernels on is equal to one. 


Theorem 3.4. Fix A G M. 


(1) Under (1.9) and (1.11), there are two nonnegative constants C and T depending on a, 
(3, 7 , and n, such that, for all (t, x) G x 

/C(t, x;\)< — Qa,fi{t, x) ^1 + F exp T t j j , (3.9) 

where a is defined in (2.6); 

(2) Under (1.9), (1.11) and the first two cases in (1.20), there are two nonnegative constants 
C and T depending on a, (3, 7 , and u, such that, for all {t, x) G x M'^, 

/C(t, x]X)>C Qa,i3{t-i x) exp ^Ai^ T t j . (3.10) 


10 









Proof. This theorem is due to Propositions 5.10, 5.11 and [4, Proposition 5.2]. 
The last set of results are the weak intermittency. 


□ 


Theorem 3.5 (Weak intermittency). Suppose that (1.9) holds and the initial data satisfy 

(3.1). 

(1) If p satisfies (2.1), then for some finite constant C > 0, 

_ 2 _ ^ 

fhp < C ^2(/3+7)-l-d/3/a pl+2(/3+7)-l-d/J/a ^ Jgp all P > 2 CVCn. 

(2) Suppose that the initial data are uniformly bounded from below, i.e., p{dx) = f{x)dx 

and f{x) > c > 0 for all x G M'^. If p satisfies ( 2 . 2 ) with |c| + | ? | 7 ^ 0, then under (1.9), 

( 1 . 11 ) and the first two cases in ( 1 . 20 ), there is some finite constant C" > 0 such that 

nip > C ^ > 2 , 

Proof. By (3.3), (3.9) and (5.11), 

I \u{t, x)\\l< Cl + C t-^ (e ^ 2 Q 2 j ^ p^t^y 

Then increase the power by a factor p/2. As for the lower bound, it holds that 

fit, a;)I Ip > I \u{t, x)\\l > + C +c^) exp f l/p" , 




thanks to (3.4) and (3.10). 

4 Fundamental solutions 

Theorem 4.1. For a G (0,2], fi G (0,2) and 7 > 0, the solution to 
+ |(-A)"/^) u{t, x) = If [fit, x)], f > 0, X G M'^, 

= Uk{x), 0 < k < \ff\ — 1, X G 


dk 


t=o 


IS 


where 


uit,x) = Joit,x) + / ds / dy fis,y) ^ Zit - s,x - y), 

Jo JR'i 


Joit,x):=Y^ / uiy]_i_kiy)d^Zit,x - y)dy 
k=o 


□ 


(4.1) 


(4.2) 


(4.3) 
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is the solution to the homogeneous equation and 




{t,x) := 


—d 


172,1 
-" 2,3 


X 




(i,i),(r/3i,/3) \ 

(d/2,a/2),{l,l),{l,a/2) i ’ 


(4.4) 


'^a,j3,'y,diti x') . tD 




^a,0,diti x) 


2,3 y 2 < 


X 




( 1 , 1 ), (/ 3 + 7 ,/ 3 ) 

(d/2,«/2), (1,1), (1 ,q/2) 


(4.5) 


and, if (d E ( 1 , 2), 


Z,g,dit,x) := ^Z„,^,rf(t,x) = TT ‘^H. 


2,1 

2,3 


X 




( 1 , 1 ), (i,« 

(rf/2,a/2),(l,l),(l,a/2) 


Moreover, 

rYcA,At,-m = 

= E,{-2-ur\(n, ah e (1.2), 


(4.6) 


(4.7) 

(4.8) 

(4.9) 


This theorem is proved in Section 4.2. For convenience, we will use the following notation 
Y(t,x) ■- Y^,„„At,x) = (4.10) 

Z*(*.i),= Z;mXi) = |z(*.i). if/3e(1.2). (4.11) 

A direct consequence of expression (4.5) is the following scaling property 

Y{t, x) = t/3+7-l-<l/3/«y r . (4.12) 


Remark 4.2. By choosing a = 2, d = 1 and (3 arbitrarily close to 2, one can see that the 
hrst condition in (1.10) suggests the condition 7 > —1. However, when 7 G (—1,0), one 
needs to specify another initial condition, namely, f{t,x)\^_^. For example (Example 
4.1 in [29, p.l38]), the differential equation 


tD^^'^gif) + g{t) = Q, (t> 0 ); 




c, 


is solved by g{t) = C — e^erfc{y/i)j. This initial condition is obscure when the driving 

term / becomes the multiplicative noise p{u{t,x))W{t,x). Hence, throughout this paper, we 
assume 7 > 0 . 


This following lemma gives the asymptotics of fundamental solutions Za^g^di^, x), x), 

and Z^ p^^{l,x) at x = 0 by choosing suitable values for 7 : 

{ \(d] in case of Z, 

/9 + 7 in case of Y, 

1 in case of Z*, when (3 G (1, 2). 
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Lemma 4.3. Suppose a G (0,2], /? e (0,2), rj G 
g{x) = (^a 

Then as x ^ 0+, the followings hold: 

' r{d-a )/2 a—d i /^( ^mm{ 2 a—d, 0 )\ 

r(r;-/3)r(o/2) i 

~r(r]-i3)r{i+d/2) log3^ + 0(l) 

2 r{l-d/a)r{d/a) . a-d\ 

ar{'n-dl3/a)r{d/2) ' > 

2r{dM + o(a;2) 


(1.1), iwh) 

(d/2,a/2), (1,1), (l,a/2) 


and d E N. Let 

X > 0. 


g{x) = < 


aV{d/2.) 


*/1? 7^ and d > a : 
*/1? 7^ and d = a : 
'll V and d < a : 

if P = V = ^- 


— ^r(l^)r(ff if (3 = g ^ 1 and d/a > 2 ■. 

r{-p)T(i+d/ 2 ) log ^ + 0{x°^) if /3 = g^l and d/a = 2-. 

2 T(l-d/a)V{d/a) . 2a-d\ 

a r(/3(l-d/a))r(d/2) ) 

0(x“) 


/3 


r(i+d/2) 

2 r(l-d/a)r(d/a) 


if ft = g ^ 1 and d/a E (1,2) : 

if (3 = g ^ 1 and d/a = 1 : 

if /3 = g 7 ^ 3 and d/a < 1 : 


Case 1, 
Case 2, 
Case 3, 
Case 4, 
Case 5, 
Case 6, 
Case 7, 
Case 8, 
Case 9, 


2 L(L-a/a)L(a/a) I 

a r(/3(i-d/Q))r(d/2) J 

where all the coefficients of the leading terms are finite and nonvanishing. 

The calculations in the proof of this lemma is quite lengthy. We postpone it to Appendix 

A.l. 


Remark 4.4. Since Dalang’s condition (1.9) implies d < 2a, the cases 5 and 6 are void 
under (1.9). Combining the rest seven cases in Lemma 4.3, we have that 


hniy'„_^,.^,d(l,a:) 


"+CXD if 7 > 0 and a < d < 2a : Cases 1-2, 
Cl if 7 > 0 and a > d = 1 : Case 3, 

C 2 if 7 = 0, /3 = 1 and a ^ d : Case 4, 

if 7 = 0, /3 7^: 1 and d < 2a : Cases 7-9, 


and 


liin Zo_;3,d(l,x) 


and when (3 G (1, 2), 


{ +CXD ii j3 ^ 1 and a < d < 2a : Cases 1-2, 
C4 ii (3 ^ 1 and a > d = 1 : Case 3, 

C2 ii (3 = 1 and a ^ d : Case 4, 


lim Z. 

x^O 


a,l3,d 


(l,x) 


+CXD if a < d < 2 q; : 
C5 if a > d = 1 : 


Cases 1-2, 
Case 3, 


(4.13) 


(4.14) 


(4.15) 


where the constants C* G M \ {0}, i = 1,..., 5, only depend on a, (3, 7 and d. Combining 
all these cases, we see that under (1.11), F(l,x) is bounded at x = 0, and under (1.12), all 
functions Z(l,x), Z*{l,x) and F(l,x) are bounded at x = 0. 
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Lemma 4.5. has the following asymptotic property as |a;| —)■ oo; 




if a = 2, 


where the nonnegative constants are 

" r(2/3+7)r(-a/2) f-J 

~ , /3{d+2-2(/3+7)) 2(/3+7)-(d+2) 

7 r-'^/2(2 - 2(2^« (2z/) 2 ( 2 -/ 3 ) if a = 2, 


and 


a = ^— ^ 6= (2-/3)/35^(2i/)?^, and c= ^ 


2-/3 


2-/3 


(4.16) 


(4.17) 


(4.18) 


Moreover, the asymptotic properties for Z{l,x) and Z*{l,x) are the same as that for Y{l,x) 
except that the argument 7 in both (4.17) and (4.18) should he replaced by \(5] — (3 and 1, 
respectively. 

These asymptotics are obtained from [20, Sections 1.5 and 1.7]. We leave the details for 
interested readers. 


Theorem 4.6. Suppose that a G (0,2], /3 G (0,2), and 7 > 0. The functions Z{t,x) : = 
Za,y,d{t,x), Y{t,x) := Ya^p^^^d{t,x) and Z*{t,x) := ZY^^^{t,x), defined in Theorem f.l, 
satisfy the following properties: 

(1) For all d ^ N and f3 G (0,1), both functions Z and Y are nonnegative. When (3 = 1, Z 
is nonnegative, and Y is nonnegative if either 7 = 0 or 7 > 1; 

(2) All functions Z, Z* and Y are nonnegative if d < 3 and 1 < /3 < a < 2. When 
l<(3 = a<2,Y is nonnegative if j > (d + 3)/2 — /3; 

(3) When d > 4, Ya^/d^Q^dif, x) assumes both positive and negative values for all a G (0, 2] and 
/3g(1,2). 

This theorem is proved in Section 4.3. It generalizes the resnlts by Mainardi et al [24] from 
one-space dimension to higher space-dimension. Moreover, in [24] only Z when f3 G (0,1] 
and Z* when (3 G (1, 2) are stndied. When /3 G (1, 2), it also generalizes the resnlts by Pskhn 
[28] from a = 2 and 7 = 0 to general a G (0, 2] and 7 > —1. 

4.1 Some special cases 

In this part, we list some special cases. 

Example 4.7. When 7 = 0 or 7 = [/3] — /3, the expressions for Z, Y and Z* in Theorem 
4.1 recover the resnlts in [9]. 


14 







Example 4.8. When a = 2, by [20, Property 2.2], we see that 




X 


and 


\2ut^ 

and, when (3 G (1,2), 


(r/31,/3) 

(d/2,1). (1,1); ’ 


W+i^P) 

(d/2,1),(1,1) y ’ 


= 7r-‘f^\x\-‘H^S ijA 


(m 

(d/2,1), (1,1) 


(4.19) 


(4.20) 


(4.21) 


In particular, for (3 G (0,1) and 7 = 0, the expressions for Z and Y recover those in [22, 17]. 
For Z 2 ^i 3 ^d, see also [21, Chapter 6]. When [3 G (1, 2), 7 = 0 and 12 = 2, the expression for Y 
recovers the result in [28]. 


Example 4.9. When a = 2 and d = 1, using Lemma A.2 and (A. 8), we see that 


z2Ai((,x) = ^ 


2x^ 


and 


>"2,/3,7,i(C 2:) = kl 
and, when (3 G (1,2), 


i\3m 

( 1 . 2 ) 


(/3+7,/3) 

( 1 . 2 ) 


^r/31-l-/3/2 


-Ml 


X 


/3/2,r/31 


^/3/2+7-1 


-M, 


4/2,/3+7 


-1 rrl.O f 2x2 




{1./3) 


t-3/2 


rM, 




X 


\fAJ 2 t3l'^ 


X 


( 1 - 2 ); \AAj2t3i‘^ I ’ 


where is the two-parameter Mainardi functions (see [4]) of order A G [0,1), 

(_l)n;^n 


Mx,,(z) := 


r(p- (n + l)A) 


, for p and z G C . 


(4.22) 


(4.23) 


(4.24) 


(4.26) 


For example, Mi/ 2 ,i{z) = exp {—z"^/A). The one-parameter Mainardi functions Mx{z) are 
used by Mainardi, et al in [24, 23]. 

Example 4.10. In [24], the fundamental solutions Za,i 3 ^d{t,x) for [3 G (0,1] and x) 

for ft G (1,2] have been studied for all a G (0,2) and d = 1. From the Mellin-Barnes 
integral representation (6.6) of [24], we can see that the reduced Green function of [24] can 
be expressed using the Fox H-function: 


^a,/3 


(x) 


H'2.1 

III ” 


X 


(1.1) , (l,/3), (l,ii=«)\ 

(1.1) ,(l,a),(l,^); ’ 


X G M, 


(4.26) 
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where a and f3 have the same meaning as this paper and 6 is the skewness: l^l < min(a, 2—a). 
For the symmetric a-stable case, i.e., 6 = 0, this expression can be simplihed using Lemma 
A.2. Hence, 


KAx) = ((kl/2) 

Therefore, their fundamental solution [24, (1.3)] 


(1/2,a/2), (1,1), (1,a/2) 


X G 


(4.27) 




2,1 




\X\ 


2^tP 


(l.l),(l,/3) 

(l/2,a/2),(l,l),(l,a/2) 


corresponds, in the case when z/ = 2, to our A when (3 G (0,1] and Z* ^ x) when 

/3g(1,2). 

Here we draw some graphs^ of these Green functions Y{t,x): see Figures 1 and 2. As f3 
approaches 2, the graphs of Y(t,x) become closer to the wave kernel \l\x\<vt/2- 



Figure 1: Some graphs of the function y2,/3,o,i(l; 3;) with n = 2, and (3 = 15/8, 5/3, 3/2, 1, 
3/4, 1/2, and 1/8 from top to bottom. 


4.2 Proof of Theorem 4.1 

Proof of Theorem f.l. Equations (4.4)-(4.9) have been proved in [9] when 7 = 0. Let / and 
'g denote the Fourier transform in the space variable and the Laplace transform in the time 
variable, respectively. Apply the Fourier transform to (4.1) to obtain 


d^uit^ 0+20 = A fit, 0 


dk 


= Ukii) , 


t=0 


^ G 

0 < A: < [,5] - 1, ^ G . 


^ The graphs are produced by truncating the infinite sum in (4.25) by the first 24 terms. In Figure 2, 
due to the bad approximations for small t when truncating the infinite sum, the graphs are produced for t 
staying away from 0. 
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(a) P = 6/5. (b) P = 3/2. (c) P = 15/8. 

Figure 2: Graphs of the Green functions x) for 1 < // < 2 for 1 < f < 6 and |a;| < 5. 


Apply the Laplace transform on the Gaputo derivative using [16, Theorem 7.1]: 

r/31-i 

C [d^ u{t, ^)] (s) = m(s, 0 ~ “fc(0- 

k=0 

On the other hand, it is known that (see, e.g., [30, (7.14)]), 

Re(7) > 0. 


£17 


Thus, 


m(s,o = + ^ i^r' 


Notice that (see [29, (1.80)]) 


01-1 _ 

UkiO + s"^/(s, 0 

k=0 


a-p 


C [t^ ^E„,^(-Af")] (s) = for Re(s) > 


Hence, 


01-1 




k=0 


from which (4.7)-(4.9) are proved. The expressions for Z and Z* in (4.4) and (4.6), respec¬ 
tively, are proved in [9]. By the fact that (see [29, (1.82)]) 

tDl {tf^-^E^,p{Xr)) = 7 e M. 

Recall that iiA“ is the Riemann-Liouville fractional derivative of order a G M (see (2.8)). 
Hence, we see that 

Y (f, x) = tD^+Z{t, x), with 6 := \/3] — // — 7, 

which can be evaluated using [20, Theorem 2.8] in the same way as in [9] for the case 7 = 0. 
This completes the proof of Theorem 4.1. □ 
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4.3 Nonnegativity of the fundamental solutions (proof of Theorem 
4.6) 

We first prove some lemmas. 

Lemma 4.11. The following Fox H-functions are nonnegative: 


(1) for all 9 G (0,1), 


H 


1,1 

2,2 


X 


(0,1), (O,0)\ 
(0,1), (0,0)7 


1 

-777; -- 77777 > 0, for X > 0 

771 + 2x^/® cos(7r6*) + x^/® 


(4.28) 


(2) for all ^ > 0 and 0 < 9 < min(l,p), 


9 X !-)■ ( |x 


( 1 , 1 ) 


> 0. 


(3) for all d gN and a G (0, 2], 


3 X H-)■ HI’I ( |x 


( 1 , 1 ) 

(d/2,o/2), (l,a/2) 


> 0 . 


(4.29) 


(4.30) 


Proof. (2) and (3) are covered by Lemma 4.5 and Theorem 3.3 of [9], respectively. As for 
(1), expression (4.28) can be found in [24, (4.38)] for the neutral-fractional diffusions. For 
completeness, we give a proof here. Because the parameters A and 6, dehned in (A.2) and 
(A.5), of this Fox H-function are equal to 0 and 1, respectively. Theorem 1.3 implies that for 
X G (0,1), 


H. 


1,1 

2,2 


X 


(0,1), (O,0)\ 
(0,1), (0,0)7 


■sp (-1)^ F(fc + 1) fc/g 
^ k\ T{-k9)T{l + k9) 

OO 

^-(-l)^sin(7rA;0)x^/^ 

k=0 

\k=o / 

1 x^/® 

77 1 + 2x^/® COs(776') + ’ 


(4.31) 


where we have applied [27, (5.5.3)] in (4.31). Similarly, when x > 1, Theorem 1.4 of [21] 
implies that 


^1,1 f (0,1), (o,0)\ ^ 'A (-1)^ _ r(fc +1) _™-(i+fc)/0 

2-2 V (0,1), (0,0)7 ^ k\ F(-(l + A;)0)F(l + (l + fc)0) 

CXD 

= ^(-l)^sin(7r(l + A;)0)x-(^+^)/® 

k=0 
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Im ( 


k^kOi^-kje 


,k=0 


= —Im 


1 + 

1 


TT 1 + 2x^/^ cos(7r6*) + 


Finally, the existence of this Fox H-function at x = 1 is not covered by Theorem 1.1 of [20] 
because A = 0 and /i = 0 (see (A.4) for the dehnition of the parameter p). In fact, as one 
can see that the series in (4.31) with x = 1 diverges. Nevertheless, we may dehne that 


^ 2,2 (1 


( 0 , 1 ), ( 0 , 0 )^ 1 
:= hm — - 


X 


1/0 


(0,1), (0,0)/ x-s-i TT 1 + 2 x^A cos(7r6') + x^/^ 27r 1 + cos(7r6*) 

This completes the proof of Lemma 4.11. 

Lemma 4.12. For /i G (0, 2], 0 G (0, 2] and d > 1, the function 


> 0 . 




(u,s) 

(d/2,1), (1,1) 


) , X > 0, 


has the following properties: 

(a) ^fd,t,,eix) = - 2 xfd+ 2 ,i,,eix). 


(b) fd,Mx) = -/= 


dz 


TT 




X^ 


-Jd+l,^,,e{z) forx > 0. 


(c) fd,^l,e{x) > 0 for all x > 0 if 9 < 2min(l,/i) and d <3. 

Proof, (a) Apply [20, Property 2.8] with k = 1, w = —d, c = 1 and a = 2 to get 


^fd,^,e{x) = X 




(0,2), (m,0) 

( 0 / 2 , 1 ), ( 1 , 1 ), ( 1 + 0 , 2 ) 


By the recurrence relation of the Gamma function, we see that 

^r(r/_~2y ^ ^ ^ d/2)F(d/2 + s) = -2F(1 + d/2 + s). 

By the dehnition of the Fox H-function, the above expression can be simplihed as 


■^fd,^^, 9 {x) = -2x (x‘^ 


(ufi) 
(( 0 + 2 )/ 2 , 1 ),( 1 , 1 ) 


) = -2x/rf+2,^,0(x). 


□ 


(b) By the dehnition of the Fox H-function, 

hh±J/Z!±£S|jA,,-2.d,, tor any 7 >-1. (4.32) 
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where the contour Lj^oo is dehned in Dehnition A.l. Assuming that we can switch the 
integrals, which can be made rigorous by writing / in the series form and applying Fubini’s 
theorem, we see that 


dz 


yj — 




-J . 


2 'ni 


ds 


r(((i +1)/2 + s)r(i + s) 

r(i-/i-0s) 


dz 


y—2s—d 


yj 


By change of variable {z/xY — 1 = y and Euler’s Beta integral (see, e.g., [27, 5.12.3 on 
p.l42]), we see that 


y—2s—d 


d;2 


X 


-2s-d /•oo 


1 2s+d 


yj z'^ — 




y 2 (l + y) 2 2 dy = 


X 


— 2s—d 


y/7iT{d/2 + s) 


r(((i + l)/2 + s) 


Note that the above integral is convergent provided that Re(2s + d) > 0, which is satished 
by choosing, e.g., 7 = Re(s) = 0 in (4.32). Therefore, 


dz 


^ -d 1 


\/ z^ — x^ 


ffd+i,ti,e{z) — X 


2 'xi 


r(d/2 + g)r(l + s) 2.^ eix) 

Til-y-Os) 2 


(c) By the recurrence in (b), we only need to prove the case d = 3. Apply Lemma A.2, 
Properties 2.4 and 2.5 in [21] to obtain 





ttI.O 

ttI.O 


7t1,0 




(m, 0/2)\ 

( 2 , 1 ) J 



(/i-0/2,6»/2) 

( 1 , 1 ) 


Then (c) is proved by an application of part (2) of Lemma 4.11. 


□ 


Proof of Theorem 4 . 6 . By comparing the Fox H-functions in (4.4), (4.5), and (4.6), We only 
need to consider the following Fox H-function: 


gix) = Hlj 


X 


(1.1), (^,/3) \ 

(d/2,a/2), (1,1), (l,a/2)y ’ 


X > 0. 


The parameter rj takes the following values 


V 


{ [/d] in case of Z, 

/d + 7 in case of Y, 

1 in case of Z*. 


(1) If /d = 1 and 7 = 0, then Z = Y and by Property 2.2 of [21], 


9 A) = 771,2 


( 1 , 1 ) 

{dl 2 ,al 2 ), (1,0/2) 


X > 0 
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which is positive by part (3) of Lemma 4.11. If /S < 1, then we can apply Theorem A.5 to 
obtain that 


9{x) = 


H 


1,1 

1,2 


( 1 , 1 ) 

(d/ 2 ,a/ 2 ), (l,a/ 2 ) 


1,0 


^1,’r 


X 


(1,1)/ t ■ 


(4.33) 


In fact, conditions (A. 12) are satisfied because 

Ai = 0, Bi = d/a, A 2 = 1, B 2 = 00. 

Moreover, = 1 and fd G (0,1) implies that = 1 — /? > 0. Hence, condition (1) of 
Theorem A.5 is satished. This proves (4.33). If /? = 1 and 7 > 0, then = A2 = 0. In view 
of condition (3) of Theorem A. 5, relation (4.33) is still true if Re(/r2) > —1 with /i2 = 1 — 7, 
i.e., 7 > 1. The two Fox H-functions in (4.33) are nonnegative by parts (2) and (3) of Lemma 
4.11. 

(2) In this case, we have that d < 3. When a = 2, by Property 2.2 of [21] and Lemma 4.12, 

9{x) = ^fy2 (1,1)) = >0. ^ > 0) 

because (d < 2 < 2min(l,?7). If a 7^: 2, by Property 2.2 of [21], we see that 


9{x) = 

As in the previous case, by Theorem A.5, we see that 


(1,1), ( 7 ?,/ 3 ), (l,a/2) 
(d/ 2 ,a/ 2 ),(l,a/ 2 ),(l,l),(l,a/ 2 ) 


HAit 


(d/ 2 ,a/ 2 ), (l,a/ 2 ) 


^2,2 


(1.1) ,(l,a/ 2 )\ ^ ^ 3^1 

(1.1) ,(l,a/ 2 )y t ~ 


X 


X > 0, 


( 1 , 1 ), (l,a/ 2 ) 

(d/ 2 ,a/ 2 ),(l,a/ 2 ),(l,l),(l,a/ 2 )y 

(4.34) 


Note that condition (A. 12) is satisfied because in this case, 

Ai = 00 , Bi = min((i, 2)/a, A2 = 1, B 2 = 0. 


When a < fd, then 


Oi = a 


/9>0, 


a*2 = 2 


a > 0, 


and condition (1) in Theorem A.5 is satisfied. When 1 < fd = a < 2, then 
02 = 2 — a > 0, a]] = Ai = 0, Re(/ii) = 1 + -—rj — 


Hence, in view of condition (2) of Theorem A. 5, the integral (4.34) is still true if 1 + d/2 — 
rj — 1/2 < —1, i.e., 7 > {d + 3)/2 — fd. 

Now, by Property 2.4 of [21], the first Fox H-function in (4.34) is equal to 


2 

a 


0 - 2,0 

-“1,2 



(»?,2 ^/q) \ 

(d/2,i),(i,i)y 


a 


21 










By Lemma 4.12 (c), we see that under the condition that — < 2min(l,?7), the hrst Fox H- 
function in (4.34) is nonnegative. This condition is satished if 1 < /5 < a < 2. By Property 
2.3 in [20], the second Fox H-function in (4.34) is equal to 

rrl,l f ^ {0,l),(0,a/2)\ 

2.2 (0,1), (0,a/2) j • 

Thanks to Lemma 4.11 (1), this function is strictly positive for t/x 7^ 0 when a G (0,2). 


(3) Now we consider the case when d > 4. The case a = 2 is covered by Lemma 25 of [28]. 
In the following, we assume that a G (0,2). By the scaling property, we may only consider 
the case t = 1. Hence, it suffices to study the following function 


g{x) = (x» 


( 1 . 1 ). (/ 3 ,/ 3 ) \ 

(d/2,a/2), (1,1), (l,a/2)y ’ 


a: > 0. 


Because a* = 2 — (3 > 0, we can apply Theorem 1.7 of [21] to obtain that 


gix) 


F((d + q;)/2) 
F(2/?)F(-«/2) 




as a: —)■ 00. 


The condition a G (0, 2) implies that F(—a/2) < 0. Thus, the coefficient of is positive. 

Hence, g can assume positive values. 

Now we consider the behavior of g{x) around zero. Because /3 > 1 and 2a < 4 < d, we 
can apply the case 6 of Lemma 4.3: 


9{x) 


F((d-2a)/2) , 

F(-/?)F(a) 




), 


as X —)■ 0+. 


The coefficient of is negative because F(—/?) > 0 for /3 G (1,2). Therefore, g{x) can 
assume negative values. This completes the proof of Theorem 4.6. □ 


5 Proofs of Theorems 3.1 and 3.2 

The proofs of Theorems 3.1 and 3.2 will follow the same arguments as the proof of [7, 
Theorem 1.2], which requires some lemmas and propositions. 

5.1 Dalang’s condition 

Lemma 5.1. Suppose that 6 > 1/2 and jd G (0,2). The following statements are true: 

(a) There is some nonnegative constant such that for allt>0 and A > 0, 

rt x2e-l 
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(b) If P < min(l,6'), then for some nonnegative constant 

ft f2e-i 

/ E‘i n{ — \w^)(\.w > Co n - , Ofl ; 

for all t > 0 and A > 0. 

Remark 5.2. When 9 = (3 = 1, then Ei{x) = e^ and thus = (2A)“^(1 — 

and (5.1) is clear for this case. 

Proof of Lemma 5.1. (a) In this case, by the asymptotic property of the Mittag-Leffier func¬ 
tion (see [29, Theorem 1.6]), for some nonnegative constants Q’s, 




w 


2(0-1) 


dtc 


= Cnt 


20-1 


(1 + \PIwfh 

-1 .., 2 ( 0 - 1 ) 


U 


(1 + XPHufh 


du 


= Cs t^^-^ 2i^i(2/3 ,29 - 1, 29- -t\^/^) 


= Ca t^^-^H. 


1,2 

2,2 


(1-2/3,1), (2(1-0),1) 

( 0 , 1 ), ( 1 - 20 , 1 ) 


(5.2) 


(5.3) 

(5.4) 


where in (5.3) we have applied [27, 15.6.1] under the condition that 9 > 1/2, and (5.4) is 
due to [20, (2.9.15)]. Notice that A = 0 for the above Fox H-function, which allows us to 
apply Theorems 1.7 and 1.11 of [20]. In particular, by [20, Theorem 1.11], we know that 


tjl,2 

“2,2 


X 


(l-2/3,l), (2(1-0),1) 

( 0 , 1 ), ( 1 - 20 , 1 ) 


0 ( 1 ), 


as x —)■ 0. 


When 0 7^ /3, by [20, Theorem 1.7], 

V (0,1), (1-20,1) /V /’ 

In particular, when 0 = /3, by Property 2.2 and (2.9.5) of [20], 


I_rl,2 

“ 2,2 


X 


(l-2/3,l), (2(1-0),1)N^ _ ^i i 

( 0 , 1 ), ( 1 - 20 , 1 ) 


) = PI I //"’■“’) = r(29 - 1)(1 + x) 


1-20 


(b) When [3 < min(l,6'), by (2.10), we know that Ej 3 ^g{—\x\) is nonnegative, hence, for 
another nonnegative constant C", one can reverse the inequality (5.2). This completes the 
proof of Lemma 5.1. □ 


Lemma 5.3 (Dalang’s condition). Let Y{t,x) = 

7 > 0. The following two conditions are eguivalent: 


a 


(i) d < 2Q(-f — min(27 —1, 0) 
P 



{t, x) with a G (0, 2], G (0, 2), and 


dyY{s, yY < oo, for all t > 0. 


23 










Proof, (i)^(ii): Fix an arbitrary t > 0. By the Plancherel theorem and (4.8), we only need 
to prove that 

fdsf <+00. (6.6) 

Jo Js.‘‘ 

Notice that d > 0 and Condition (1) together imply that /3 + 7 > 1/2. Thus, we can integrate 
ds first using Lemma 5.1 (a). Then it reduces to prove that 

X _|_ |^| 2 «+| min(27-l,0)^ + 00 , 



which is guaranteed by (i). 

(ii)=^(i): The case (3 G (0,1] can be proved in the same way as above by an application of 
Lemma 5.1 (b). The case (3 G (1,2) is trickier. Fix t > 0. Denote the integral in (5.5) by 
I{t). Then by change of variables, 

pt poo 

I{t) = C 

Jo Jo 

Note that the double integral is decoupled. The integrability of ds at zero implies that 
2(/9 + 7) “ 1 ~ > 0, which is equivalent to 


d<2a + |(27-l). (5.6) 

y 

By the asymptotics of ?/) at +cx) (see, e.g.. Theorem 1.3 in [29]; note that the 

condition (3 G (1,2) is used here), we see that there exist |/o > 0 and some constant C > 0 
such that 

C 

Ely+^i-y) > for all y > y^. 

Hence, the integrability of dy at zero and infinity implies the following conditions: 

d/a>d and d/a —3<1. (5.7) 

Combining (5.6) and (5.7) gives (i). This completes the proof of Lemma 5.3. □ 

5.2 Some continuity results on Y 

This part contains some continuity results on Y. All the results proved in this part will be 
used in the proof of Theorem 3.2. In particular. Proposition 5.4 will be used to prove the 
Holder continuity (Theorem 3.3). 

Proposition 5.4. Suppose a G (0,2], (3 G (0,2), 7 > 0, and (1.9) holds. Then Y{t,x) = 
ya,i 3 ,'y,dit,x) Satisfies the following two properties: 
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(i) For all 0 < 9 < {Q — d) A 2 and T > 0, there is some nonnegative constant C = 
C{a, f3, 7, u, 6, T, d) such that for all t G (0, T] and x,y E 


drdz (Y {t — r, X — z) — Y {t — r, y — z))'^ < C 


X 



(5.8) 


(a) If (3 <1 and 7 < \(Y\ — [3, then there is some nonnegative constant C = C{a,(3,'y, iy,d) 
such that for all s,t E with s <t, and x E 

[ dr f dz (Y {t — r, X — z) — Y {s — r, X — z))'^ < C{t — s)2(/3+7)-i-<i/3/a ^ 

Jo Jr^ 

and 


dr / dzY^t-r,x-z) < C{t - s)2(/3+7)-i-rf/3/a 


(5.10) 


Proof, (i) Fix t > 0. By Plancherel’s theorem and (4.8), the left hand side of (5.8) is eqnal 
to 


{2t^Y 

1 

1 

< 


dr (f - / dfE\p^^{-2-^n{t-rY\fY) 

d^ (1 - cos(^ • {x - y))) [ dr{t- rY^h+z-d (-2"V(f - r)^|^|") 

^ Jo 


(21,) 


d( (1 - cosK ■ (a: - »))) 


X _|_ |^|2«+f min(0,27-l) 


where we have applied Lemma 5.1 in the last step (see also the proof of Lemma 5.3). Denote 
0 := 2a + I min(0, 27 — 1). Becanse 1 — cos(x) < 2 A {x^/2) for all x G M, we only need to 
bonnd 


^^ 2A[|x^ <C^| \x-y\-'^ 


.V2 


u 


d+l 


1 + iei 


(1 + |x — y\~^u)^ 


du 


+ \x-y\ 


e-d 


IV 2 


e+i-d 


du 


u 


The second integral on the right hand side of the above ineqnality is finite provided that 
Q > d, which is Dalang’s condition. By [27, 15.6.1], for some constant C > 0, 


-72 


——- 1 .g dn = C 2aFi(0, 2 + d, 3 + d; -a/2|x - y\~Yi 

(1 + |x — y\~'-u)^ 


which is trne nnder the condition that d + 3>d + 2>0. By [20, 2.9.15], 


-72 


(1 + \x 


n 1 ti~\~ 1 / 

-= C"Hl'^ (V2\x - y\-^ 

X — y\ ’ V 


(-1-71), (1-0,1) 

(0,1), (-2-d,l) 
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Since A = 0, by [20, Theorem 1.7], for all 9 G (O,min(0,2 + d)), 


-x/2 


u 


d-\-l 


/ nil I 1 = 0{\x-y\ ), as|a:-?/|^0. 

Jo [I + \x - y\-^u)'^ 

Combining these cases, we have proved (i). 

(ii) Denote the left hand side of (5.9) by I. Apply Plancherel’s theorem and use (2.10), 

I = crdr f de 

Jo JM.d. 

-{s - (- 2 -V(s - r)^|e|“) 

Then by Lemma 5.5 below, 

I = CCi [ dr [(t - r) 2 C+T'-i)-'^A« + _ ^) 2 (/ 3 + 7 -i)-d/ 3 /aj _ 2 C f dr[{t-r){s-r)f+^-^H{r), 

Jo Jo 

where 

H(r)^ [ Bft^+,(- 2 -V(«-rf| 5 r)£;^,,j+,(- 2 -V(s-r)'>|J|‘>)dJ, 

JR'* 

By (2.10) and t > s, 

//(>■)> / Bj^+,(- 2 -V((-rf|^r)d{ 

JR'' 

= (t _ ^)-2C+7-1) / (i _ 

= (t — f Y{t — r,yydy 

JR'^ 

= Q(t-r)-"^/“, 

where in the last step we have applied Lemma 5.5. Because /? + 7 < 1, we see that 
[ dr[{t — r){s — H{r) > ( dr(f — 

Jo Jo 

Denote p := 2[(3 + 7 ) — 1 — d(3/a. Note that p > 0 is implied by Dalang’s condition (1.9). 
Therefore, 

/ < [tP _ (t - s)P + sP- 2{{t - sf - sO] < ^(t - s)C 


This proves (5.9). As for (5.10), by a similar reasoning, we have 


ca 


/ dr dzY^{t-r,x-z)< CC, / dr (f - = - Ht - s)p . 

Js Jr^ Js p 

This completes the proof of Proposition 5.4. 


□ 
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The following lemma is a slight extension of [26, Lemma 1] from the case where 7 = 1 —/? 
to a general 7 . 

Lemma 5.5. Assume that d < 2a, (3 G (0,2), and 7 > 0. Then 

[ x)Ax = C, 


for all t > 0, where 


a: = 


r(d/2)(27rz/)'^/2 




(5.11) 


Proof. By Plancherel’s theorem and (4.8), 

^ 2 (/ 3 + 7 - 1 ) 


Y{t,xydx = 




( 2 vr)'^ jK<i 
2T^d/2 ^ 2 (/ 3 + 7 - 1 ) 


r(d/ 2 ) ( 2 vr)^ 7 o 

2T^d/2 ^2(/3+7-l)-<i/3/a /-oo 

f777 (aTp I 


r‘^-^E|_^+.^(-2-VfV)dr 




Note that by the asymptotic property of the Mittag-Leffler function ([29, Theorem 1.7]), the 
last integral is finite if d < 2 q;. □ 

The corresponding results to the next Proposition for the SHE, the SFHE, and the SWE 
can be found in [7, Proposition 5.3], [8, Proposition 4.7], and [6, Lemma 3.2], respectively. 
We need some notation: for r > 0, a > 0 and (t, x) G x denote 

Bt,x,T,a '■= {{t', x') G X : 0 < f' < f + r, \x — x'\ < a^ . 

Proposition 5.6. Suppose that (3 G (0,2) and 7 G [0, [/d]—/d]. Then for all {t,x) G M;]_xM'^, 
there exists a constant H > 0 such that for all (t',x') G and all s G [0,P) and ?/ G 

with \y\ > A, we have that Y (t' — s,x' — y) < Y (t + 1 — s, x — y). 


Proof. Without loss of generality, assume that u = 2. 

Case I. We hrst prove the case where a = 2. The proof here simplihes the arguments 
of [4, Proposition 6.1]. Fix {t,x) G x M'^. By the scaling property and the asymptotic 
property of E, we have that 

Y{t + 1 — S,X — y) / t' — S I2; _ y|a / ]j\^x _ y\c \ 

Y{t' — s,x' — y) + 1 — s/ \x'— y\°- ^\(d —s)^E 2 x _ 5^/372 J ’ 


as \y\ —)■ 00 , where the constants a, b and c are dehned in (4.18), and 


m 



- /3-iy + 


aft 

Y' 
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Notice that 


t + 1 — s t + 1 — t' t + 1 — t + 1 1 

-= IH-> IH-> -- = IH- 

t' — s t' — s t' t + 1/2 2 t + l 


If /(/9) < 0 , then 

If /(/?) > 0, then 



/t + 1 — 
V t'- s 



\fm 


> 1 . 


(5.12) 


t' — S \ ■f ^ ft' — s 
t + 1 — s J yt + l 


\fm 


(i + l)-l/OTIexp(|/(/3)| log(i'-s)). 


The rest arguments are the same as the proof of [4, Proposition 6.1]. We will not repeat 
here. 


Case II. Now we consider the case when a G (0,2). By the scaling property and the 
asymptotic property of Y, we have that 


Y{t + 1- s,x -y) _ / t'-s \ 12/3^/ 

Y{t' — s,x' — y) \t + l — s) \\x — y\) 


as \y\ —)■ oo. Because /3 > 1/2 and 7 > 0, we see that 1 — 2/3 — 7 < 0. Hence, by (5.12), 


t' — s 
t + 1 — s 


1-213-u 


ft + 1 — s 

V t'-s 


2 / 3 + 1/-1 

> 



1 

2t + 1 


213+u-l 

> 1 . 


On the other hand, 

\x' - y\ ^ \y\ - |x^| ^ \y\ - (\x' - a:| + | 3 :|) ^ |y| - 1 + |x 
\x-y\ ~ \x\ + \y\ ~ \x\ + \y\ ~ \x\ + \y\ 


as \y\ —)■ 00 . Therefore, we can choose a large constant A, such that for all \y\ > A and all 
{t',x') e 1 / 2,1 and s G [ 0 ,T], 


Y{t + l-s,x-y) ^ ^ 

Y (f — s,x' — y) 

This completes the proof of Proposition 5.6. □ 

Proposition 5.7. For all {t, x) G M+ x 1 < /3 < 2 and 7 G [0, 2 — /3], we have 

lim / / dsdi/ {Y {t' — s,x' — y) — Y {t — s,x — y))"^ = 0. 

Proof. This proposition is a consequence of Proposition 5.6. The proof follows the same 
arguments as the proof of [4, Proposition 6.4]. □ 
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5.3 Estimations of the kernel function 1C 

Let G : M+ x i—)■ M with d G N, h > 1 be a Borel measurable function. 

Assumption 5.8. The function G : M+ x i—)■ M has the following properties: 

(1) There is a nonnegative function Q{t,x), called reference kernel function, and constants 
Co > 0, cr < 1 such that 

G(Mr<^e(i.n. tor all (M) e R+ X K". (5.13) 

(2) The reference kernel function Q{t,x) satishes the following sub-semigroup property, for 
some constant Ci > 0, 

/ dy Q {t,x — y)Q (s, y) < Gi Q {t s, x) , for all t,s > 0 and x G (5-14) 


Assumption 5.9. The same as Assumption 5.8 except that the two “<” in (5.13) and 
(5.14) are replaced by “>”. We call the property (5.14) with “<” replaced by “>” the 
super-semigroup property. 

Proposition 5.10. Under conditions (1.9) and (1.11), the function Y{t,x) satisfies As¬ 
sumption 5.8 with the reference kernel Ga,/ 3 {t,x) defined in (3.7), two nonnegative constants 
Go and Gi, depending on (a,/9,7, z/, d), and a defined in (2.6). 

Proof. The proof is similar to that of [4, Proposition 5.8]. We hrst note that a < 1 is implied 
by Dalang’s condition (1.9); see (2.7). 

Case I. We hrst consider the case where a = 2. In this case, 

Notice that 

^^^>L/3J+1, for/3 G (0,1) U (1,2), 

and when (3 = 1, the constant b dehned in (4.18) reduces to l/(2z/), which is bigger than 
l/(4z/). Hence, by (4.16) and (4.13), we see that 

Y{t,xY Y{l,yy ^ ^ 

sup — -- = sup —— -- =: Co < oo- 

{t,x)(£R+xR‘i t "^02,git, X) y^Rd y2,g{l,y) 

Note that in the application of (4.13) in the above equations we have used the fact that 
Dalang’s condition (1.9) implies d < 2a. When /3 = 1, we see that 

^2,i(t,a;) = (4z/7rf)“'^/%xp , 
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and hence, Ci = 1 and (5.14) becomes eqnality in this case. When (3 G (0,1), 


= (4z/7rt^) j < JJ (4z/7rt^) ^^%xp 


Z=1 


\Xi 


4z/t^/2 / 


Then by Lemma 5.10 of [4], for some nonnegative constant Ci, 

/ ^2,/3(^ -s,x- |/)^2,/3 (s, i/)dy < Cl Q2,/3(t, a:). 

When /3 G (1,2), 

^2,/3(t,a;) = (4z/7rt^)"'^^%xp = ^2,i(t^,a;). 

Hence, by the semigronp property for the heat kernel, 

[ Q 2 ,fiit - s,x- y)G 2 ,fsis, y)dy = ^2,i((t - sf + /, x) < a;), 

jRd 

where in the last step we have applied the ineqnalities: 

<{t-sf+ s^ <t^, for/3 g (1 ,2). 

Hence, in this case, Ci = Therefore, Assnmption 5.8 is satished. 

Case II. We now consider the case where a ^ 2. In this case, 

c„, 


Ga,y (j^i a^) 




= Gp(f^/“,x), 


where Gp{t, x) is the Poisson kernel (see [32, Theorem 1.14]). By the scaling property and the 
asymptotic property of Y{l,x) at 0 and oo shown in (4.13) and (4.16), for some nonnegative 
constant C, 

(J -j^y+'y—l-dy/a fj -j.y+'i-l-dp/a 

Y{t,x)< 


< 


(1 + t-2/3/«|x|2)('^+“^/^ (1 + t-2/3/«|x|2)(^+^)/^’ 

where the second ineqnality is dne to {d + a)/2 > (d + l)/4, which is eqnivalent to d> l — 2a. 
Hence, 


sup 


Y(t,xf Y(l,yf ^ ^ 

— t;— T, —\ ^0 < OO- 


(t,x)eR+xR‘^ t ^Qa,pidiX') j/eRd ^q:,/3(1) ?/) 

Then, it is ready to see that 

Qas{t,x) = Gp{t^^°‘,x) 

By the semigroup property of the Poisson kernel, we have that 


GaAt -s,x- y)Ga,/3(s, y)dy = Gp (s^/“ + {t - s)^/", x). 


(5.15) 
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Then use the inequalities 


^/3/a ^ _ ^^/3/a < 2^-^/“^^/“ if /?/« < 1, 

2l-/3/a^/3/a <g/3/a ^ _ g)/3/a < ^/3/a 


(5.16) 


to conclude that for some constant Ci > 0, 



GaAt -S,X - y)GaA^^ VAy ^ = ClGaAt, A- 


This completes the proof of Proposition 5.10. 


□ 


Proposition 5.11. Under (1.9) and the first two cases of (1.20), the function Y{t,x) sat¬ 
isfies Assumption 5.9 with the reference kernel Qa,i 3 it,x) defined in (3.8), two nonnegative 
constants Co and Ci, depending on (a,/3, 7 , z/, d), and a defined in (2.6). 

Proof. The proof is similar to the proof of the previous proposition. We have several cases. 

Case I: When a = 2 and fi G (0,1), by the scaling property (4.12), and the asymptotics at 
zero and inhnity in (4.13) and (4.16), we see that 


t ‘^^2,/3(t,a;) _ ^2,/3(1,2/) _ 1 


(5.17) 


By the semigroup property of the heat kernel. 



where the last inequality is due to 

A <{t- sf P for/3e (0,1). 

Case II: When a < 2 and /? G (0,1 V a), by (4.12) and (4.16), 

CQaAt^x) QaA^^y) 1 ^ , 

"""P . Yit xP ^ yri vV 

By the super-semigroup property can be proved in the same way from (5.15) and (5.16). □ 

5.4 A lemma on the initial data 

Lemma 5.12. For all compact sets K C x Mfi, 


sup ([1 + Jq] * (^5 A < 

{t,x)£K 


under the following two cases: 
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(1) Both (1.9) and (1.11) are satisfied and the initial data satisfy (3.1); 

(2) Both (1.9) and (1.12) are satisfied and the initial data belong to Ma,/s (I^)- 
Proof. In both cases, the kernel function }C{t,x) has the following upper bound 

IC{t, x; A) < CiGaAt^ A ■ 


Part (1) is clear because 


([1 + J2] * /C) (t, x) < (1 + a) (1 * /C) (f, x) 


Ciil + Ct) 



where cr < 1 (see (2.7)). The proof of part (2) requires more work. The case when a 
proved in Lemma 6.7 of [4]. The proof for a G (0, 2) is similar to that of Lemma 4.9 
Let 


G«,/3,d(t,x) = TT ^|x| 


2,1 

2,3 


X 


2 “ ^vth 


( 1 , 1 ), (»?,/ 3 ) 

(d/2,a/2), (1,1), (l,a/2) 


= 2 is 
in [8]. 


with 7] = [/d] in case of Z and r) = 1 in case of Z*. Hence, we need only consider jQ(t,x) = 
(|/i| * G{t, •))(2^)- By the asymptotic properties both at inhnity and at zero (see Lemma 4.5 
and Remark 4.4), we have that for all t G (0,T], 


G{t,x) = f^-^-'^^/“G'(l,f-^/“x) < 


(J ^r\-\-dg/a (J ^ri-l-d/S/a^^ y rpyP 


1 + 


< 


1 + lx 


\d-\-<y. 


Thus, for s G (0, f], 

Jfis,y) < V tffi 

where 

The rest of the proof follows line-by-line the proof of part (2) of Lemma 4.9 in [8]. This 
completes the proof of Lemma 5.12. □ 


5.5 Proof of Theorem 3.2 

Proof of Theorem 3.2. The proof is the same as the proof of [4, Theorem 3.1], which in turn 
follows the same six steps as those in the proof of [7, Theorem 2.4] with some minor changes: 

The proof relies on estimates on the kernel function /C(f, x), which is given by Proposition 
5.10. 

In the Picard iteration scheme, one needs to check the L^(fI)-continuity of the stochastic 
integral. This will guarantee that the integrand in the next step is again in V 2 , via [7, Propo¬ 
sition 3.4]. Here, the statement of [7, Proposition 3.4] is still true by replacing in its proof 
[7, Proposition 3.5] by either Proposition 5.4 for the slow diffusion equations or Proposition 
5.7 for the fast diffusion equations, and replacing [7, Proposition 5.3] by Proposition 5.6. 
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In the first step of the Picard iteration scheme, the following property, which determines 
the set of the admissible initial data, needs to be verihed: for all compact sets K C M_,_ x 

sup ([l + Jq] * /C) (t, x) < +00. 

For the SHE, this property is proved in [7, Lemma 3.9]. Here, Lemma 5.12 gives the desired 
result with minimal requirements on the initial data. This property, together with the 
calculation of the upper bound on }C{t,x) in Theorem 3.4, guarantees that all the Lp{Q)- 
moments of u{t, x) are hnite. This property is also used to establish uniform convergence of 
the Picard iteration scheme, hence LP(r2)-continuity of (t,x) i—)■ I{t,x). 

The proof of (3.4) is identical to that of the corresponding property in [7, Theorem 2.4]. 
This completes the proof of Theorem 3.2. □ 


5.6 Proof of Theorem 3.1 


Proof of Theorem 3.1. The proof of Theorem 3.1 is similar to that for Theorem 3.2. Because 


Ct = sup 

(s,a:)G[0,t] 


|Jo(s,a:)| < cx), for all t > 0, 


the Picard iterations in the proof of Theorem 2.4 [7] give the following the moment formula 


\u{t,x)\\l < 2.J{t,xf + 


+2C'f 


1-k JCp] (f, x). 


Note that the function (1 x JCp) {t, x) is a function of t only. For convenience, we denote it as 

H{t-,X):=[ ds [ dylC{s,y,X). (5.18) 

Jo JRd. 

Therefore, we need only to prove that H{t-, X) is finite, which is proved in Lemma 5.13 below. 
This completes the proof of Theorem 3.1. □ 

Lemma 5.13. For all a G (0,2], jj G (0,2), 7 > 0, and d G N, under Dalang’s condition 
(1.9), we have that 

H{t; X) < exp (^CX^^ t^ , 

for allt>0 and A G M, where a is defined in (2.6) and C is some constant depending on a, 
(3, 7 and d. 


Proof. By Lemma 5.5, 

(1 X £o)(t, X) < Q £ ds s2(/3+7-l)-<i/3/a ^ 

where (Pji is defined in (5.11) and 9 = 1 — a. Note that 6* > 0 is guaranteed by Dalang’s 
condition (1.9). Now we claim that, for n > 0, 


(lx£„)(f,a;) < 


C', 


n+lp ^(n+l)0 

F((n +1)0 + 1) 


(5.19) 
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of which the case n = 0 is just proved. Assume that (5.19) holds for n. By the above 
calculations, we see that 


(1 * Cn+i) (t, x) < ^ I ds (f - ^ 


r((n + 1 ) 6 * + 1 ) Jq 


T{{n +2)6 + 1) 


Therefore, 




n=0 


Then apply the asymptotic property of the Mittag-LefHer function (see, e.g., [29, Theorem 
1.3]). This completes the proof of Lemma 5.13. □ 


A Appendix: Some properties of the Fox H-functions 

In this section, we follow the notation of [20]. 

Definition A.l. Let m, n,p, g be integers such that 0<m<g,0<n<p. Let 0^,6* G C 
be complex numbers and let aj, /3j be positive numbers, i = 1, 2, • • • ,p; j = 1, 2, • • • ,q. Let 
the set of poles of the gamma functions r{bj + /3js) doesn’t intersect with that of the gamma 
functions r(l — Oj — Ojs), namely. 


—bj — ^ , 1 f 1 — Ui T ^ , I 

bji = —^-, / = 0, 1, • ■ ■ f n 1 ^ ^-’ k = 0,l,--- } = 




(y,i 


for alH = 1,2, • • ■ ,p and j = 1, 2, • • ■ ,q. Denote 


mn( ^ _ nr=i + «is) nr=i r(i ais) 

■ nr=.+i r(a, + a,.) nu+i r(i - ’ 


The Fox H-function 


TTm,n/\ — Tjm^n 
^p,q [Z) = 


is defined by the following integral 


PQ 


1 

27ri 


(ai,ai) 

(&i,A) ■■■ 

(®P) ^p) 

i.bg,Pq) _ 

j H;‘^{s)z-‘ds, 

z eC, 


(A.l) 


where an empty product in (A.l) means 1, and L in (A.l) is the infinite contour which 
separates all the points bji to the left and all the points Oik to the right of L. Moreover, L 
has one of the following forms: 

(1) L = L_oo is a left loop situated in a horizontal strip starting at point — cxd + and 
terminating at point —oo + for some —oo < (j)i < (j )2 < oo 

(2) L = L+oo is a right loop situated in a horizontal strip starting at point +cxd + and 
terminating at point oo + icj )2 for some — oo < 01 < 02 < oo 
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( 3 ) L = Lj-yoo is a contour starting at point 7 — ioo and terminating at point 7 + zoo for 
some 7 G (—00, 00) 


According to [20, Theorem 1.1], the integral (A.l) exists, for example, when 

q p 

A := (3j — Q!i > 0 and L = L_oo, 

j=i i=i 


(A.2) 


or when 


a* := Qfj — Qfj + I3j — (3j > 0 and L = Li^^. 

i=l 2=72+1 j=l j=m-\-l 


(A.3) 


The following two parameters of the Fox H-functions (A.l) will be used in this paper: 

q p 


/i = ^ ^ Oi + 

j = l i=l 


P-Q 
2 ’ 


and 


(A.4) 


i=l j=l 


j ■ 


Lemma A.2. For 6 G C and (3 > 0, there holds the relation 


Tj2+m,n 
^P,2+q ' ^ 


(6,/3),(b+l/2,/3),(b,-,/3Pi,, 


ol-26 rz Tjl+m,n 

= 2 Vtt 


and 


Tjm.n 

^p,2+q I ^ 


{ct-l ± 

(b,/3),(6+1/2,/3) 


) = 4-+-I/2 ^l+rn^n ^ 


(26,2/3),(6j ,/3j )i,q- 

(ttj l,p 

(26,2/3) 


Proof. Let 


. n^i r(fei ++s) nr=i r(i - a* - a^s) 
n{s) = 


ULn +1 r(a* + ttis) n?=™+i r(i - - + s) 

By the dehnition of the Fox H-function, 


lj2+m,n 
^p,2+q I ^ 


(A.5) 


[ - / 'H{s)T{b + (ds)T(b + 1/2 + (ds)z ®ds. 

(6,/3),(6+l/2,/3),(6,-,/3Pi,,; 2 Tri J^ 


By the duplication rule of the Gamma function [27, 5.5.5 on p. 138] 

F(z)F(z + 1/2) = +^2^-2^F(2 z), 2z ^ 0, -1, -2,..., 

we have that 

1 


(A.6) 


lj2+m,n 
^p,2+q > ^ 


{o-i ■jO'i) 1 

(6,/3),(6+1/2,/3),(6, 


) / W(s)r(2() + 2/)s)(4+)-*ds. 

>pj)T-,q/ 2m r 


Then apply the dehnition of the Fox H-function. The second relation can be proved similarly. 

□ 
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Here are some direct consequences of this lemma: 


H 


l+m,n 

l+p,l+g 


(26,2/3), (6j,/3j)i,q/ 


o26—1^—1/2 Tjl+ni,n 

2. vr ^p,i+q 



(l/2+6,/3),(6,,/3,)i.J ’ 


Trl+m^n f + 

1+P,1+'? V (26,2/3), (6,,/3,)i,, y 


<-)26—1„—1/2 Tjl+m,n 

2 tt 



(<^i,Q^i)l,p \ 

(6,/3),(6,,/3,)i,J ’ 


H 


m,l+n 

l+p,l+q 


(1/2+6,/3), (ai,ai)i,p\ 
(bj ,0j)i^q, (26,2/3) y 


4^7r 


1/2 Tjm,n 

^p^+q 



(ai,aj)l^p \ 

(6,,/3,)i,„(6,/3)j ’ 


hr 


m,l+n 

l+p,l+(j 


(1/2+6,/3), (ai,ai)i,p\ 
(bj ,0j)i^q, (26,2/3) y 


4^7r 


1/2 Tjm,n 





Remark A.3. In [4], the Green function Gp{t,x), which corresponds to x), is 

represented using the two-parameter Mainardi function of order A G [0,1) (see (4.25)). By the 
series expansion of the Fox H-function ([20, Theorem 1.3], which requires that A = 1 —A > 0), 
one can see that 


[z I j/")) , A e |0,1), 

By Property 2.4 of [20], the above relation can also be written as 

= Ag[ 0,1). 


(A.7) 


(A.8) 


Remark A.4. Another commonly used special function in this setting, such as in [28], is 
Wright’s function [36, 35, 37] (see also [23, Appendix F]): 


OO ^ 

:= ^ for A >-1,/I e C. (A.9) 

^ k\T{fi + Xk) 


We adopt the notation 0 that is used by E. M. Wright in his original papers. By (2.9.29) 
and Property 2.5 of [20], 


0(A,/x;z) 



HoS [z 

(i)i)i (1+-^—Mi-^) 

A” (‘ 

(p-A,-A)\ 

(1,1) J 


if A > 0, 
if A G ( — 1; 0] . 


Comparing (A.7) and (A. 10), we see that 


Ma,^( 2 :) = 0(-A,/i - A;z), for A G (0,1]. 


(A.IO) 


(All) 


The following theorem is a simplihed version of Theorems 2.9 and 2.10 in [21], which is 
sufficient for our use in the proof of Theorem 4.6. 
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Theorem A.5. Let (a^, Ai,/ii) and (a^, A2,/i2) be the constants (a*, A,/i) defined in (A.3), 
(A.2) and (A.4) for the following two Fox H-functions: 


I ^ 


, tjM,N ( 

, , and rip'r, [x 


(ci-7i)i,Q 


respectively. Denote 


Ai = min 

l<i<n 


Ao = min 


1 - Re{ai) D _ 


(Xi 

Re{cj) 


Bi = min 

i< ■' 

B 2 = min 


l<j<m ’ 

1 — Re{di) 

l<i<N 


l<j<M ' ■ l<i<N Si 

with the convention that min( 0 ) = +00. If either of the following four conditions holds 

(1) > 0 and > 0; 

(2) = Ai = 0, Re{pii) < —1 and > 0; 

(3) 02 = A2 = 0, Re{fi2) < —1 and al > 0; 

(4) = Al = 0, i?e(/xi) < —1 and = A2 = 0, Re{fi 2 ) < —I, 

and if 


Al + Bi > 0, A 2 + B 2 > 0, Al + A 2 > 0, Bi + B 2 > 0, 


(A.12) 


then, for all z > 0, x G M, 


-^p+P,(?+Q 



{cti^<i^i)n-\-l,p 


f 


n r 




X 

t 


Proof. By Property 2.3 of [21], 


{di,Si)i^p\ dt 
{ci,7j)l,Q/ t 


H 


M,N 

P,Q 



{di,5i)i^P 

(ci>7i)l,Q 


ttN,M 

^Q,P 


X 


(l-Cj,7j)l,Q 

R-di,Si)l,P 


(A.13) 


If condition (1) holds, one can apply Theorem 2.9 of [21] with p = 0, a = 1, w = 1/x, 
and with the following replacements: N —>■ M, M^N,P^Q,Q^P, Cj^l — Cj, 
dj —)■ 1 — dj. If either of conditions (2)-(4) holds, we apply Theorem 2.10 of [21] in the same 
way. Note that the parameters p for both Fox H-functions in (A.13) are equal. □ 


A.l Proof of Lemma 4.3 

Proof of Lemma 4 ■ 3. Let 


f{x) = Hl’l (^x 


[ipAwh) 

(d/2,a/2), (1,1), (l,a/2) 
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Then g{x) = x ^f{x°'). Let 


H, 


d,a,l3,r] 


,s : = 


T{d /2 + as/2)T{l + s)T{-s) 
r{r] + f3s)r{-as/2) 

Denote the poles of r(l + s) and T{d/2 + as/2) by 

21+ d 


A := { —(1 + /c) : A; = 0,1, 2, ■ ■ ■ } and B ; = 


a 




respectively. According to the dehnition of Fox H-function, to calculate the asymptotic at 
zero, we need to calculate the residue of Hd^a,/ 3 ,ri{,s)z~^ at the rightmost poles in A U i?. 
Because a* = 2 — /3 > 0, all the nigh cases are covered by either (1.8.1) or (1.8.2) of [21]. 
The notation h*i below follows from (1.3.5) of [21]. 


Case 1. Assume that f] ^ /3 and d/a > 1. In this case, the rightmost residue in A U i? is 
at s = — 1 and it is a simple pole. Hence, 


r((<i-a)/2) 
r(, - / 3 )r(a/ 2 ) 


> 0 , 


and 

f(x) = h* 2 QX + O ^ Qg ^ ^ g_|__ 


Case 2. Assume that tj ^ f3 and d/a = 1. The rightmost residue in A U H is at s = —1 
and it is of order two. Hence, 

Res{Hd 4 ,f},n{s)x~") = lim [{s + l)‘^Hd 4 ,^^rj{s)x~^]' 

= ym ([(s + - (s + l)^ffj,j,^,,(s)logx) x~‘ 

where we have used the fact that r(a:) has simple poles at x = —n, n = 0, 1 ,..., with residue 
. Therefore, 


Case 3. Assume that rj ^ (3 and d/a < 1. The rightmost residue in A U H is at s = —d/a 
and it is a simple pole. Hence, 


^ 2 V{1 — d/a)V{d/a) 

^ aV{r]-dp/a)V{d/2) 


(A.14) 


where the nonnegativity is due to the fact that rj > (3 > (3d/a. Therefore, 

f{x) = + 0 (a;““(A+ 2 )/«.i)) = + 0(x), as X ^ 0_ 
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Case 4. Assume that rj 


/3 = 1. By Property 2.2 of [21], 


f{x) = hI’I 


( 1 . 1 ) 

(d/2,«/2), (l,a/2) 


Hence, 


and 


_ 2T{d/a) 
aT{d/2) 


7 ^ 0 , 


f{x) = h\^x^l^ + as a: ^ 0+. 


Case 5. Assume that rj = (3 ^ 1 and d/a > 2. The rightmost residue in AUB is at s = —1, 
but this residue is vanishing because hms_)._i l/r(/d + /3s) = 0. The rightmost nonvanishing 
residue in A U H is at s = — 2 and it is a simple pole. Hence, 


h 


* 

21 


r((d- 2 a)/ 2 ) 

r(-/3)r(a) 


Hence, 

f{x) = hl^x^ + 0(a;™°(3,d/A)^ as X ^ 0+. 


Case 6. Assume that rj = fd 1 and d/a = 2. As in Case 6, the rightmost nonvanishing 
residue in A U H is at s = — 2, and it is of order two. Then 


Res(hfrf,rf/2,/9,/3(s)x ^) = ^hmj(s + 2)^Trd,d/2,/3,/3(s)x *]' 

= lim ([(s + 2fHd4/2,i3A^)]' -i^ + 2)^hfd,d/2,/3,/3(s) logx) x~ 


s=-2 


s —>--2 

= Cx^ + 


T{-(3)T{l + d/2) 


logx. 


Therefore, 




Case 7. Assume that r] = (3 ^ 1 and d/a & (1)2). As in Case 6, because h* 2 Q = 0, the 
rightmost nonvanishing residue in A U i? is at s = —d/a, and it is a simple pole. Hence, 

/(x) = h^Qx'^^^ + O(x^), as X —)■ 0+, 

where is defined in (A. 14) with rj replaced by f3. 

Case 8. Assume that r] = /3 ^ 1 and d/a = 1. The rightmost nonvanishing residue in 
A U H is at s = —1, and it is of order two. Hence, 

Res{Hd4,0^l3{s)x~^) = lim [(s + l)^i7rf,d,/3,/3(s)a:"^]' 

S = —1 —1 

= lim [Hi( s)''H 2(s) + Hi(5)7^2(s)'- 7^1 (5)7^2(s) logx]x“®. 
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where 


•Hi(s) = (s + l)2r((l + s)rf/2)r(l + s) and 'H 2 {s) = 
As calculated in the proof of Lemma 7.1 of [9], we have that 


r(-s) 


r(/3 +/3s)r(-ds/2) ■ 


2 ("l + ci2 2 

'Hi(-l) = lim 'H\{s) = - = lim -- -r = 

s—^—i d s—1 -)- s)ci/2)(l As) d 

H 2 {- 1 ) = lim ni{s) = 0, 

S-S--1 




= lim 


r(-s) 


s=-l 


-ir(-ds/2) Vr(/3(i + s)) 

r(i) ^im + s)) 


r(d/2).^-i r(/3(i + s)) 

/9 

r(d/ 2 ) ’ 

where 'ip{z) is the digamma function and the last limit is due to (5.7.6) and (5.7.1) of [27]. 
Therefore, 

= r(i + rf/2) 


and 


fix) = 




r(i + d/ 2 ) 


x + 0{x^), as X —)■ 0+. 


Case 9. Assume that rj = /3 ^ 1 and d/a < 1. The hrst nonvanishing residue in A n i? is 
at s = —d/a and it is a simple pole. Hence, 

f{x) = hlQx‘^^°‘ + (9(x), as X —)■ 0+, 

where is dehned in (A. 14) with rj replaced by (3. This completes the whole proof of 
Lemma 4.3. □ 
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